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AN ADAPTIVE CONTROLLER FOR SPACE BASE 


SPIN PLANE ORIENTATION 

By Willard W. Anderson and Ralph W. Will 
Langley Research Center 

SUMMARY 

An adaptive attitude control system for the precise earth orientation of an artificial- 
gravity space base is presented. The adaptive logic forms a performance index based on 
the control torques required to hold a constant orientation with respect to the orbit plane 
in the presence of environmental torques. The adaptive control system acts to minimize 
this performance index and thus acquires an orientation with respect to the orbit plane 
where the environmental disturbance bias torques precess the space base at the orbit 
regression rate. This precession permits the space base to operate at a fixed orientation 
with respect to the orbit plane with no requirement for steady-state attitude-control fuel. 

It is also shown that this condition is stable and that any angular offset from this point 
results in a stable coning motion of the space base spin axis with respect to the orbit 
plane. Simulation results are presented for both the adaptive control system perfor- 
mance and the space base motions. 


INTRODUCTION 

Manned orbiting space base concepts, such as that shown in figure 1, involve large 
rotating elements for artificial -gravity modules and zero-gravity hub modules for engi- 
neering and scientific experiments. Pointing experiments such as astronomical tele- 
scopes and earth resource packages require that the zero-gravity hub and thus the 
rotating element spin axis be stabilized. Arbitrary pointing and stabilization of the 
spacecraft spin or momentum vector require significant amounts of reaction jet fuel or 
prohibitively large momentum storage devices because of the extremely large angular 
momentum associated with the space base rotating element. A method of stabilizing the 
space base by utilizing the external gravity gradient and aerodynamic disturbances would 
eliminate the need for continuous attitude control. 

One means of achieving such space base passive stabilization is to tilt the rotating 
element momentum vector out of the orbit plane so that the bias components of the exter- 
nal disturbances would precess the spacecraft at the orbit regression rate and therefore 
maintain a constant orientation with respect to the orbit plane. Selection of the proper 



spacecraft orientation for this operation would be accomplished by a relatively simple 
adaptive controller. The initial acquisition using adaptive control techniques would be 
made by using reaction jets to maneuver the entire space base. After acquisition, base 
control would be possible by use of a control moment gyro (CMG) system mounted in the 
nonrotating element. The authors wish to acknowledge parallel efforts of K. L. Lindsay 
of Manned Spacecraft Center. 

These considerations were used to synthesize adaptive control logic, an adaptive 
system performance index, and a space base control law to acquire the proper orientation 
for passive control of the space base. 


SYMBOLS 


A matrix of direction cosines 

a arbitrary maneuver magnitude 

ajj matrix elements 

Ci,C2,Cs quasi-linear constants 
D maximum drag force 

E single -rotation matrix 

F aerodynamic drag force 

G unit bias torque 

H angular momentum 

I principal inertia 

i orbit inclination 

J gravity gradient torque 

K adaptive controller gain 

k local-vertical direction cosines 
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aerodynamic torque 
longitude 

orbit precession rate 
space base spin rate 
center-of -pressure coordinate 
body torque 
time 

unit vector 

coordinate axes 

state variable 

stopping condition bound 

equatorial-ecliptic inclination 

complex coning vector 

performance index 

orbit period 

time constant 

Euler angle set 

low-frequency Euler angle set 

base-spin-vector offset angle 


base-spin-vector cone angle 



fi 


orbit angular rate 


spin frequency 

S2 0 linear base precession frequency 

w base body rate 

Subscripts: 


a 

actual data 

b 

body 

c 

control 

d 

aerodynamic drag 

h 

heliocentric 

i 

inertial 

o 

orbit 

s 

sun 

V 

local vertical 

x,y,z 

axes of space base 

Conventions: 


matrix 

column matrix or vector 

absolute value 

vector gradient operator 


[ ] 
( ) 
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A change in a variable 

( ) average value for one base revolution 

« )> average value for one base orbit 

Bars over symbols denote vectors. Dots over symbols denote derivatives with 
respect to time. Primed symbols denote spatial derivatives. 

SPACE BASE AND CONTROL SYSTEM CONFIGURATION 

Figure 1 shows a possible space base configuration to house 50 to 100 men and a 
scientific laboratory in the 1980 time period. The space base shown consists of two 
basic elements: a rotating artificial-gravity living quarters section ("disk") and a zero- 
gravity hub which contains all pointing experiments and the logistics system docking 
facilities. The two elements are connected by a low-frequency spring system contained 
in the compliant hub module which isolates the hub and its critical pointing experiments 
from the disturbances due to wobbling and imbalance of the rotating section. The spring 
mechanism also allows the hub to be maneuvered with respect to the rotating element 
spin axis and is sized so that the natural frequencies of the two rigid-body systems lie 
above the gravity gradient cyclic frequency (twice orbit frequency) but below the rota- 
tional frequency. Therefore, the bodies behave as a single rigid body from the cyclic 
gravity gradient point of view but the hub is isolated from the higher frequency "disk" 
rotation inputs. 

The zero-gravity hub is stabilized and pointed by a control-moment-gyro (CMG) 
system. This orientation method includes rolling the entire hub about the base spin axis 
to provide a second degree of freedom for the solar cell array and solar experiments 
module. The pointing experiments are mounted on individual gimbaled platforms which 
provide the precise stabilization involved and the pointing capability with respect to the 
hub required to work with various targets such as the earth's surface, the sun, and stellar 
phenomena. A large space shuttle vehicle (12 -man capacity) could be docked to the end of 
the hub module and would represent the most radical change in configuration inertia, and 
associated gravity gradient disturbances experienced by the base during normal operation. 

The rotating section of the space base is controlled and reoriented by a reaction jet 
system located on the ends of the 200-foot rotating arms. The reaction jets are pulse- 
modulated to fire during parts of the spin cycle as shown in figure 2 to produce torques 
about particular axes. It should be noted that the firing regions in figure 2 are refer- 
enced to the "low-frequency body axes" corresponding to the system with respect to which 
the base is controlled. When the rotating base arms are within these regions, the jets 
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may be fired for the proper amount of time and in the proper sense to produce the desired 
precession rate for reorientation or control. 


Physical System Description 

The basic coordinate system used for inertial reference is a geocentric equatorial 
coordinate system X^Y^Zj and is shown in figure 3. The Xj-axis is arbitrarily chosen 
to pass through the vernal equinox. This axis is collinear with the X^-axis of a helio- 
centric ecliptic coordinate system Xj^Yj^Zh, the rotation about the X-axis being the 
equatorial-ecliptic inclination e (e ~ 23.5°). A sun- line unit vector u g is defined by 
the heliocentric longitude l s . The space base orbit plane Xq,Y 0 is also shown; it is 
defined by the equatorial longitude l and the orbit inclination i. Orbit precession l 
is taken as -4.5° per day for an inclination of 55° and a circular orbit of 250 nautical 
miles (463 km). (See ref. 1.) Space base attitude, relative to the orbit coordinate sys- 
tem, is defined by the Euler angle set 4> 0 ,S 0 ,^ 0 for the 1,2,3 rotation sequence in fig- 
ure 4. Space base position in orbit is given by the orbit plane longitude l\y Also illus- 
trated in figure 4 are a local vertical unit vector u v and a free-stream air velocity unit 
vector u^. 

Space base attitudes relative to orbit coordinates are required for control purposes 
but are not an inertial set, since the orbit regresses. A second set of Euler angles 
!/A,$i,0i for a 3,2,1 rotation sequence, relative to the geocentric equatorial coordinate 
system, is used to describe the physical system. These Euler angles are found by 
assuming the space base spin momentum magnitude to be constant for the range of fre- 
quencies considered and by solving the following statement of conservation of angular 
momentum: 

T b - cn b x H (1) 

where T b is the total body torque, aJ b is the body rate, and H is the constant magni- 
tude spin momentum. This formulation excludes the higher frequency dynamics associated 
with base flexibility, rotation, and nutation, but does give an accurate low-frequency char- 
acterization of an actual space base provided the spin momentum is large and remains 
alined with the Zb-axis. The body rates found by solution of equation (1) are transformed 
to Euler rates and integrated to yield the set ^,^,0.. The transformation of body rates 
to Euler rates, taken from reference 2, is 




The Euler angle set 0 o ,0 o ,i// o are found by taking the matrix of direction cosines 
defining the body in inertial coordinates [Ab{] and transforming it to a matrix of direc- 
tion cosines defining the body in orbit coordinates jA^oJ. The elements of this matrix 
can be used to solve for the three Euler angles as shown in appendix A. (See ref. 3.) 


The environmental torques considered in this report are gravity gradient torques 
and aerodynamic torques. These torques are derived in appendix B and their final form 
excludes the "high” frequency components associated with base revolution. From appen- 
dix B, the gravity-gradient torques in body coordinates are 


<J X > = -3S22 ky k z (<I x>y ) - I z ) 
<J y ) = -3G 2 k x k z (l z - <I x , y >)> 
<Jz> = 0 


(3) 


where 


{ k } “ [ A bo]{ u ^ 


The aerodynamic torques, also in body coordinates, from appendix B are 


< L x> r cp,y F z “ r cp,z F y 

(Ly) = r C p jZ F x - r C p )X F z \ 

< L z> = r C p,x F y " r cp,y F x^ 

where 

(F) . (<D>} T [A bo ](u^ 


(4) 


Of these two sets of environmental torques, only the gravity- gradient torques are 
assumed to have steady-state components (bias). These bias torques can be shown to 
vary linearly with body-to-orbit attitude, for small body-to-orbit angles, and for the 
example configuration of this report can be characterized as follows (appendix B, 
eqs. (B7)): 

«Jx» - -G4> 0 1 

«J y » * -Ge 0 I 

«J Z » - 0 ( 


G « | n 2 (l z - <I x ,y>)J 


The following linear (small 4> 0 ,® 0 >4 / 0 J analysis is included to characterize the pas- 
sive coning of the space base spin vector in terms of angle and frequency. 
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It is is assumed that \p Q is initially zero and remains so since «J Z )> « 0, 
referring to figure 5, the X^-axis and Y^-axis body rates are 


GBr 


w x H 


y H 




( 6 ) 


J 


The orbit plane regression rate components along the X^-axis and Y^-axis are zero and 
i sin(i + <p^j, respectively. Transforming these rates to find rates for the body relative 
to orbit yields 


* 

v o h 


/ . G0„ 

e Q = -l sin(i + 

or expanding (sin i + <p^J and rewriting yields 




(7) 


where 


+o* il o*o = ^o , | 

e o = g^o 


* _ Hi sin i/G 


( 8 ) 


1 + 


Oq — 


III cos i 
G 

(!cAV2 


(9) 


G 1 + 


HZCi 


H 


A complex solution to equations (8) for an arbitrary set of initial conditions (<£ o (0), 
0 o (O)), valid when Sl Q » |z cos ij or |HiCj/G| « 1, is shown in figure 6 where 


£ = 4 >* + 0**e -in o t (10) 

and where 

C = f ? o «» - + [ e o (°>] 2 } l/2 

This relationship means that, within the stated assumptions, a space base spin 
vector will cone about an axis which is stationary relative to the regressing orbit plane, 
the angle and period of coning being given by <p ** and S2 0 , respectively. This report 
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r 


describes a steepest descent adaptive controller with and without learning the type that 
would be utilized on a spinning space base to seek out and aline the base with the station- 
ary vector and, thus, maintain the base with a fixed attitude relative to its orbit plane, 
without the use of steady -state fuel. 

Control System Description 

The function of the adaptive controller (fig. 7) is to find that space base orbit orien- 
tation where no control torques (fuel) are required to maintain the orientation. The adap- 
tive controller performance index p must thus be defined as a function of control torques 
which, in turn, are functions of the orientation angles. The problem is reduced to two 
degrees of freedom, since the spinning space base considered averages out the Xb-axis 
and Yb-axis inertia difference. (See appendix B.) However, to insure a two-degree-of- 
freedom problem, a Zb-axis controller is used to maintain at zero. The perfor- 

mance index p is an even function of the average Xb-axis and Yb-axis control torques 
for one-half orbit. The half-orbit does remove the influence of cyclic gravity -gradient 
components, since their period is one-half orbit, but it does not necessarily remove the 
influence of cyclic aerodynamic torques, since their period is an entire orbit. The 
shorter time period is obviously desirable and no adverse effect from using it occurred 
during simulation. The index is defined as 

p{<P 0 ’ 9 o) =< T c,x> 2 +<Tc,y> 2 (11) 

where 



and r is the orbit period. 

The performance index p is measured by holding the base at a fixed orientation 
and computing the integrals of the torques required to hold that orientation. If the atti- 
tude control system errors are small during this period, these integrals represent an 
accurate measurement of the external surface and body torques acting on the base, the 
components of torques required to maintain orbit regression rate not being included. By 
orienting the base with a single -degree -of -freedom maneuver (for example, <p and 
remeasuring p, an estimate of the variation of p with this degree of freedom (for exam- 
ple, 9 p/9<^ 0 ) is found. A similar maneuver for the remaining degree of freedom 9 Q 
allows an estimate of dp/d0 o . The two maneuvers allow an estimate of Vp. An adap- 
tive control technique is to allow these maneuvers to have amplitudes which are propor- 
tional to the measured gradients. This method is called the method of steepest descent. 
(See ref. 4.) The successive orientations are thus 
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( 12 ) 


x new x old ~ ^ V| ° 

where 



A stopping condition is met when | Vp|^ is less than a fixed amount 6. No constraint 
boundaries are included and the gain K is held constant. 


The maneuver and hold logic used first computes desired Euler rates from the new 
Euler command angles 


^o,c ~ T^(^o,new ~ ^o) 

^o,c = ^(^o,new " ^o) ) 

*o,c = ^K) 

J 

and then converts these angles to body rates 

w x, c = 0 o ,c cos ^O cos ^o + ^o,c sin 

w y, c = <?o,c cos " ^o,c cos e o sin / 

w z, c = *o,C + <Vc sin 0 O 


( 13 ) 


( 14 ) 


and commands the base reaction control system to produce these body rates. The system 
continues this iterative motion until the stopping condition is met. 


The steepest descent logic, as outlined, operates from point to point in state vector 
space and does not retain information prior to the previous iteration. A learning pro- 
cedure was included to use prior information in order to reduce the time required to 
reach the steady-state orientation. Past maneuver amplitudes are used to form two 
Taylor series expansions from which future maneuver amplitudes (&<P Q C ,A0 O c j can be 
estimated. These two Taylor series are both solved iteratively for a point of mini- 
mum p (no commanded maneuver) near a point determined by linear extrapolation 
only. These expansions for command maneuvers are 




( 15 ) 


O = A0 o>a +A ^ 


where the primes denote approximate spatial derivatives. 


J 
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An algorithm is used to compute and test successive maneuvers for appropriate 
stability during the initial steepest descent iterations and when these maneuvers meet 
certain stability tests, cj) Q and 0 O maneuvers are commanded according to the solu- 
tions of equations (15). After acquisition of this assumed point of minimum p, the 
method of steepest descent logic is used to maneuver to within the stopping condition 
bounds of the exact point of minimum p. 

SIMULATION RESULTS 


The space base low-frequency dynamics, low-frequency environmental torques, 
coordinate transformations, and adaptive control logic were programed on a Control 
Data 6600 digital computer. Table I includes a summary of assumed inertia, spin rate, 
and angular momentum characteristics of the example space base configurations con- 
sidered during the simulation. Also shown in table I are the offset angles and frequen- 
cies predicted by the linear analysis (eqs. (9)) for the four configurations, and their 
agreement with simulation runs. 


TABLE I.- EXAMPLE CONFIGURATION SUMMARY 
[i = 55°; 1 - -4.5°/day; altitude, 250 nautical miles (463 km)j 


Case 

<Ix,y> 


z 

N, 

H 

<b 

0> de R 

f2 0 , rad/sec 


slug-ft2 

kg- 

slug-ft2 

kg- m2 

rpm 

ft-lb-sec 

N-m-sec 

Linear 

Simulation 

Linear 

Simulation 

Low speed 

193 X 10 6 

262 x 10 6 

320 X 10 6 

CD 

O 

X 

'cr 

eo 

0.2 

6.75 x 10 6 

9.15 x 10 6 

1.133 

1.15 

38 X 10~ 6 

36.2 x 10' 6 

Design 

193 

262 

320 

434 

2.0 

67.5 

91.5 

12.95 

13.0 

3.55 

3.09 

Shuttle docked 

243 

330 

320 

434 

2.0 

67.5 

91.5 

23.60 

25.1 

2.03 

*.89 

Growth version 

353 

479 

640 

868 

2.0 

135 

183 

11.28 

11.4 

4.03 

3.59 


*Dynamics are too nonlinear for agreement. 


The space base design configuration rotates at 2 revolutions per minute to produce 
0.2 5g at the center line of the living quarters module. The low speed configuration 
represents the basic space base spinning at 0.2 revolution per minute. This case was 
considered to demonstrate that the adaptive logic is not sensitive to spin speed for a 
particular inertia configuration. A docked shuttle spacecraft adds an inertia load of 
50 x 10® slug-ft^ (68 x 10® kg-m^) to the space base X- and Y-axes (Z-axis change 
neglected) and represents a radically different inertia distribution. The growth version 
of the space base consists of four rotating arms rather than two for the design configu- 
ration. This configuration would represent a total increase in base crew contingent or 
power generation capability. 

Typical simulation output is shown in figure 8, which is a time history of space 
base attitudes and environmental and control torques for the 0.2-rpm acquisition case 
of figure 9. The spacecraft attitudes <p Q and 6 Q are measured with respect to the 
orbit plane as defined in figure 4. Also shown are the gravity -gradient disturbance 
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torques which are referenced to the space base body axes. The orbit precession 
angle and the control system torques are also shown. A corresponding cross-axis 
trace of the spin axis angles, cp 0 and d Q , with respect to the orbit plane is shown 
in figure 9. This trace was generated simultaneously with figure 8. The spin axis 
started normal to the orbit plane (<p 0 ( 0) = 0 o (O) = Oj and was allowed to precess under the 
external disturbances. The motion agrees well with the linear prediction represented in 
figure 6, measured values of (p* and £2 0 agreeing to within 5 percent of the linear 
predictions of table I. The cyclic gravity-gradient disturbances (at twice orbit frequency) 
are shown on the higher speed sections of figure 8 and their influence on base motion is 
evident on the cross-axis trace. The adaptive system was then activated and started to 
acquire the point where the external disturbance bias will maintain the base orientation 
constant with respect to the precessing orbit. The control torques during this period 
are low-level cyclic torques required to counteract the external disturbances and high- 
level step torques required to perform the reorientation maneuvers commanded by the 
adaptive logic. The adaptive logic may be terminated at any point as shown in figure 9 
and the spacecraft allowed to precess with the resulting error or coning angle decrease. 
The adaptive control system was reactivated and completed the acquisition of the desired 
orientation. The system was shut off automatically when the stopping condition was met. 
The external disturbances again precess the spacecraft in a small cone with respect to 
the orbit plane, the cone angle being a direct function of the stopping condition. Super- 
imposed on this low-frequency coning motion is the precession caused by the cyclic 
gravity-gradient and aerodynamic disturbances. 

Figure 10 shows a cross-axis trace for a similar run involving the 2-rpm design 
space base configuration. This configuration, having ten times more angular momentum, 
is precessed through approximately ten times the previous cone angle by the external dis- 
turbances. The measured values of <p * and 0 (;) agreed to within 13 percent of those 
listed in table I. The only modification to the adaptive system for the 2.0-rpm configura- 
tion required was a change in the initial maneuver angles from 0.1° to 1°. This change 
is necessary because the gradient is not as large for this case and requires a large angle 
for accurate measurement. Figure 11 shows adaptive control system acquisition for the 
2-rpm case from a variety of initial conditions with respect to the orbit plane. No diffi- 
culty was experienced by the adaptive controller for these conditions. 

The docking of a shuttle spacecraft to the space base, as shown in figure 1, will 
radically change the base inertia distribution, the external disturbance profile, and the 
orientation at which the spacecraft will be precessed at orbit precession rate by the 
external disturbances. Figure 12 shows a cross-axis trace involving the 2-rpm config- 
uration where acquisition has been accomplished by the adaptive control system and a 
shuttle vehicle docking occurs which upsets the stable situation. The adaptive controller 
is left off for a period following the docking to illustrate that the external environment 
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will precess the spacecraft in a new cone. The measured value of <f>* agreed with that 
in table I; however, <f>* is too large for the linear analysis to predict S2 0 , where the 
measured value was approximately one -half that given in table I. However, when the 
adaptive system is activated, no difficulty is experienced in acquiring the new desired 
orientation with the combined configuration. 

The growth version of the space base, which has four rotating arms as shown in 
figure 1, represents another extreme configuration variation for the adaptive control sys- 
tem. Acquisition of the desired orientation from a variety of initial conditions for this 
configuration is shown in figure 13. Values of <p * and S2 0 measured agreed to within 
11 percent with those of table I. Minor controller modifications, consisting of doubling 
the maximum reaction jet control torque level (using a second set of reaction jets located 
on the second set of arms) and decreasing the adaptive gain k in proportion to the 
increase in performance index, resulted from the increased inertia differences. This 
technique permitted the controller to accomplish acquisition without difficulty. 

The simplicity, predictability, and reliability of the method of steepest adaptive con- 
troller are attractive; however, the number of orbits required to acquire the steady-state 
orientation is large and results in a large expenditure of fuel. For the design configura- 
tion acquisitions shown in figure 11, the average fuel consumption was approximately 
850 lbm (385 kg) (for an assumed specific impulse of 250 lbf-sec/lbm (2440 N-sec^/kg)) 
whereas the average number of orbits was 70. Total spin-up fuel for this configuration 
was 1120 lbm (508 kg). 

Figure 14 shows design case acquisitions using the learning procedure discussed 
previously. The acquisition sequence begins by establishing a steepest descent pattern. 
When the ratio of successive maneuvers meets a stability criterion of the control algo- 
rithm, the space base is maneuvered directly to the computed point of minimum p. The 
method of steepest descent logic is then reactivated, if required, and the base is maneu- 
vered until the stopping condition is met. For the design case acquisitions shown in fig- 
ure 14 the average number of orbits required was 50 which represents considerable 
improvement over identical initial condition runs without the learning procedure. 

CONCLUDING REMARKS 

In the preceding analysis, equations for the "low-frequency" space base dynamic 
system have been developed and the system simulated. The results of this analysis and 
the simulation rims allow certain concluding remarks concerning the behavior of the sys- 
tem, within the boundaries of all connected simplifying assumptions. 

First, the basic uncontrolled precession of the space base, with respect to the 
regressing orbit, resulting from the bias components of the gravity torques describes 
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a stable cone of predictable amplitude, offset, and frequency. Superimposed on this basic 
motion is the stable base coning resulting from the cyclic gravity gradient and aerody- 
namic torques. 

Secondly, a simple and inherently reliable method of steepest descent adaptive con- 
troller has been synthesized and simulated which will acquire that axis in orbit space 
which minimizes the passive gravity -gradient-bias coning to any desired level. The con- 
troller can be turned off when the desired orientation is reached and no further expendi- 
ture of fuel is required. 

Finally, the benefits of a simple learning technique have been demonstrated to illus- 
trate the reduction of time and fuel that is possible if the complexity of the control algo- 
rithms is not limited. 

Langley Research Center, 

National Aeronautics and Space Administration, 

Hampton, Va., June 3, 1970. 


14 


APPENDIX A 


COORDINATE SYSTEM TRANSFORMATIONS 

Space base attitude relative to a noninertial orbit coordinate system is a required 
input for the adaptive control computer. The transformations for obtaining this informa- 
tion used in the simulation are summarized in this section. The simulation transforma- 
tions are analogous to those required to transform an inertial reference unit's attitude 
output to an orbit reference for control purposes. 

Given the space base inertial attitude as a set of 3,2,1 rotation sequence Euler 
angles of magnitude i/q, 0j, and <p^, respectively, and using the notation of reference 3, 
the matrix of direction cosines required to transform a body vector to an inertial vector 
is 

A ib = E 3 (-*i) E 2 (-0 i ) 


cos cos 0^ -cos 0. sin 0| + sin 0^ sin cos 0^ 

cos sin 0^ cos 0 i cos i + sin 0. sin sin 0j 

-sin sin 0 i cos 9^ 


sin cb^ sin 0^ + cos 0^ sin 9 ^ cos 0^ 
-sin 0^ cos 0^ + cos 0^ sin 6 ^ sin 0^ 
cos 0 i cos 


(Al) 


where 


(}i = MOb 

The matrix of direction cosines required to transform an inertial vector to an orbit vector 
is 


Aoi = Ej(i) E 3 (i) 


cos l 

-cos i sin l 
sin i sin l 


sin l 0 

cos i cos l sin i 

-sin i cos l cos i 


(A2) 


where 

Oo - MO. 

The product of these two matrices Aq^ = AqjA^ yields the matrix of direction cosines 
to transform a body vector to an orbit vector. The inverse of this product, A^ 0 = A^ 
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APPENDIX A — Continued 

(or transpose since it is orthogonal) can be used to compute the orbit reference Euler 
angles shown in figure 4. This computation is accomplished by equating the matrix 
A bo(0i> 0 i>^i>M) to the matrix A bofoo^o’^o)’ This latter matrix is 

Abo = E 3(^o) E 2( 0 o) E l(^ 0 ) 

cos ip 0 cos e Q cos sin d Q sin 0 Q + sin i i/ Q cos <£ q -cos sin 6 0 cos (p Q + sin sin 

= -sin ip Q cos e 0 -sin \jj Q sin 9 Q sin <f> Q + cos i p Q cos <f> Q sin \p Q sin 8 0 cos <p Q + cos \p Q sin 4> Q (A3) 

sin e 0 -cos 8 0 sin <p Q cos 9 Q cos 

where 

Ob = iAbojOo 

The following procedure for computing d Q , and (p 0 can be verified by inspection of 

equation (A3) or by referring to reference 3: 

i/' = tan - * 

0 a ll 

5 1 = a 13 sin ^o + a 23 cos ^o 
C 1 = 3-22 cos + a l2 sin 

5 2 = a 31 

c 2 = a 33 C l " a 32 S l 

i S 1 

<jf> = tan " 1 -L 

r o c 1 

e a = tan-' 

where the elements of Ajj 0 are 

a H = cos l cos cos + sin l cos 0j sin ip. 

a 12 = -cos t s t n l cos ^i cos + cos i cos l cos 0^ sin i/a - sin i sin 0^ 
a ^3 = sin i sin l cos 0j cos - sin i cos l cos 0j sin i/a - cos i sin 0j 
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APPENDIX A - Concluded 

&21 = cos z(-cos sin 0^ + sin 0^ sin 0j cos 0j) + sin l {cos 0^ cos 0^ 

+ sin sin 0j sin ip^j 

a 22 = -cos i sin z(-cos 0^ sin 0j + sin 0^ sin 0j cos 0^ + cos i cos z(cos 0^ cos 0^ 
+ sin 0. sin 0^ sin xp^j + sin i sin 0^ cos 0j 

a 23 = s ^ n * s ^ n ^( -cos 4>i sin 0j + s i n sin cos 0j) - sin i cos ^(cos 0^ cos 0^ 

+ sin 0^ sin &i sin 0jJ + cos i(sin 0^ cos 9^j 

agj = cos z(sin 0^ sin 0^ + cos 0^ sin 0^ cos 0^ + sin z(-sin 0. cos 0^ 

+ cos 0^ sin 0^ sin 0^ 

a 32 = " cos i sin l (sin 0j sin 0^ + cos 4>j_ si - n ®i cos 0i) + cos i cos z(-sin 0^ cos 0^ 
+ cos 0^ sin 0^ sin 0jj + sin i cos 0^ cos 0j 

agg = sin i sin z(sin 0j sin 0 i + cos 0. sin 0j cos 0^ - sin i cos Z (-sin 0^ cos 0j 
+ cos 0^ sin 0 i sin 0.^ + cos i(cos 0 i cos 0^ 
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APPENDIX B 


ENVIRONMENTAL TORQUES 


The following set of gravity gradient body torques are from reference 5 and have 
been rewritten in the nomenclature of this report: 


Jjj = -3fi 2 kyk z (ly “ I^j 

Jy = -3£2 2 k x k z ^I z - I x ) \ (Bl) 

J z = -3fi 2 k x ky^I x - Iy) 


where the direction cosines k x , ky, and k z define the local vertical unit vector at the 
base center of mass in body coordinates, {k} = [AbJ]^u^ . 

Since the base rotation frequency is much greater than the orbit frequency, this con- 
tribution can be eliminated for the low-frequency cases considered as follows: 


Expanding 


yields 


where 


® - [Abo] 



k x = cos ^(Cj) + sin ^ 0 (C 2 ) 
k y = cos ^ O (o 2 ) - sin i// Q (Ci) ) 
= C3 


Cj = cos 0 O cos l\y + sin 0 o sin <p Q sin 
C 2 = cos sin 

C 3 = sin 6 0 cos Zj, - cos 6 Q sin sin 


(B2) 


The terms Cj, C 2 , and C 3 are considered to be constant for the period of one base 
rotation. Also the Z-axis attitude is considered to have a low-frequency component 

(constant for one rotation) in addition to the rotation or spin component, or 


^o = K + 


(B3) 
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APPENDIX B - Continued 


where fi b is the base spin velocity. By referring to figure 4, it is noted that is 
the final Euler rotation, and thus can be manipulated without changing (f> 0 or 9 Q , as is 
shown in figure 15. Rewriting the gravity gradient torques for the low-frequency body 
axes (x b ,Y b ,Z b ^, of figure 15 by using equations (Bl), (B2), and (B3), averaging for one 
rotation of the body, and omitting the asterisks, yields 


(J x ) — 3R 2 (c 2 cos - Cj sin ^ 0 ) (^3)( Ox,y ) - *z) 
<J y > = -3S2 2 (+C! cos \f/ Q + C 2 sin ^ 0 )(C 3 )(l z - <I x ,y>)> 
<J Z ) = 0 


(B4) 


where 


I* + Iv 
<!x,y) = g 

From inspection of equations (B2) , equations (B4) can be written as 


<J X > = -3fi 2 kyk z (<I X> y> - I z )" 
<Jy> = -3fi 2 k x k z (l z - <I X ,y ) ) > 

<J Z ) = 0 


(B5) 


By assuming to be zero and <p Q and 0 o to be small, equations (B4) sim- 
plify to 


<J X > ~ -3R 2 (0 O sin Z b cos l h - <p Q sin 2 Z b )(<I x>y > - I z ) 
<J y ) * -3 Sl 2 (e Q cos 2 l b - <p Q sin l b cos f b )( I Z - <I x ,y>)) 
<Jz) ~ 0 


(B6) 


If it is assumed also that <p Q and 0 O remain essentially constant during one-half 
orbit, then 


«J X » * -| fi2(l z . <I X> y>)0 o 
«Jy» « fi 2 (l z ~ <Ix,y>)0Q / 


(B7) 



APPENDIX B - Concluded 


Aerodynamic torques result when the center of pressure of a body does not coincide 
with the center of mass. If a vector in body coordinates r cp defines the position of the 
center of pressure relative to the center of mass, then the torque resulting from an aero- 
dynamic drag force F also written in body coordinates, will be 

L = ^cp X F 

The drag force F can be modeled as having a direction parallel to the free- stream air 
velocity and a magnitude proportional to the scalar product of a vector (D) whose 
components are the averaged maximum drag forces along the low-frequency body axes, 
and the free-stream unit vector u^, transformed to body coordinates. A matrix equation 
expressing this relation is 


(<L >} = [r C p] (<D >} T [A bo ] (u d ) 


(B8) 


Expanding equation (B8) yields 


(L x ) = r C p, y F z - r C p jZ F y 

(Ly) = r C p ;Z F x - r C p )X F z \ 

<Lz> = r cp,xFy “ r cp,yFx 


(B9) 


where 


F x 


F z 


[cos i// Q cos 0 o sin - [cos \)/ Q sin 0 o sin <p Q 
|jsin cos 9 Q sin l b + [sin xj/ Q sin 0 o sin <p Q 
[sin 9 0 sin b b + cos 9 0 sin cos Zb|(D z ) 


+ sin i// Q cos cos J b ]<D x > 
- cos \p Q cos 0 Q )cos ^b](Dy> 


It should be noted that these aerodynamic torques averaged for one orbit are zero 
(if it is assumed that the orientation is constant, r cp is constant, and that the average 
maximum drag forces D x , Dy, and D z are constant for the entire orbit), or 


«L X » = «Ly» = «L Z » = 0 
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Figure 7.- Adaptive controller block diagram. 





























Figure 15.- Gravity gradient torque components. 
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